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Abstract 

We study the differential branching ratio, forward-backward asymmetry, CP-violating asym- 
metry, CP-violating asymmetry in the forward-backward asymmetry and polarization asymme- 
tries of the final lepton in the B — » ^ + ^ _ decays in the context of a CP softly broken two Higgs 
doublet model. We analyze the dependencies of these observables on the model parameters by 
paying a special attention to the effects of neutral Higgs boson (NHB) exchanges and possible 
CP violating effects. We find that NHB effects are quite significant for the r mode. The above- 
mentioned observables seems to be promising as a testing ground for new physics beyond the SM, 
especially for the existence of the CP-violating phase in the theory. 

PACS numbers: 12.60.Fr, 13.20.He 

1 Introduction 

Although CP violation is one of the most fundamental phenomena in particle physics it is still one 
of the the least tested aspects of the the Standard Model (SM). Before the start of the B factories, 
CP violation has only been measured in the kaon system. Very recently, the observation of CP 
violation in the B-meson system have been reported by the e + e~ B factories [ 1 1 providing the the 
first test of the SM CP violation. In the near future, more experimental tests will be possible at the 
B factories and possible deviations from the SM predictions will provide important clues about 
physics beyond it. This situation makes the search for CP violation in B decays highly interesting. 

Interest in CP violation is not limited to particle physics; it plays an important role in cos- 
mology, too. One of the necessary conditions to generate the matter-antimatter asymmetry ob- 
served in the Universe is -in addition to baryon number violation and deviations from the thermal 
equilibrium- that the elementary interactions have to violate CP. In the SM the only source of CP 
violation is the complex Cabibbo-Kobayashi-Maskawa (CKM) matrix elements which appears 
too weak to drive such an asymmetry [2], giving a strong motivation to search for new physics. In 
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many cases, extensions of the SM such as the Two Higgs Doublet Model (2HDM) or the super- 
symmetric extensions of the SM are able to supply the new sources of CP violation, providing an 
opportunity to investigate the new physics by analyzing the CP violating effects. 

Being a FCNC process, B — ► X s ^ £ + £~ decays provide the most reliable testing grounds 
for the SM at the loop level and they are also sensitive to new physics. In addition, B — > 
Xd£ + £~mode. is especially important in the CKM phenomenology. In case of the b — > s£ + £~ 
decays, the matrix element receives a combination of various contributions from the intermediate 
t, c or u quarks with factors V t bV t * s ~ A 2 , V c bV* s ~ A 2 and V u bV* s ~ A 4 , respectively, where 
A = sin 6c — 0.22. Since the last factor is extremely small compared to the other two we can ne- 
glect it and this reduces the unitarity relation for the CKM factors to the form V t b V t * + V c bV* s « 0. 
Hence, the matrix element for the b — > s£ + £~ decays involve only one independent CKM factor 
so that CP violation would not show up. On the other hand, as pointed out before QJ @1, for 
b — ► d£ + l~ decay, all the CKM factors V tb V t * d , V cb V* d and V ub V* d are at the same order A 3 in 
the SM and the matrix element for these processes would have sizable interference terms, so as 
to induce a CP violating asymmetry between the decay rates of the reactions b — > d£ + £~ and 
b — ► d£ + £~ . Therefore, b — ► d£ + £~ decays seem to be suitable for establishing CP violation in B 
mesons. 

We note that the inclusive B — ► X s £ + £~ decays have been widely studied in the framework 
of the SM and its various extensions 151- 1221 . As for B — ► Xd£ + £~ modes, they were first 
considered within the SM in [3] and [4]. The general two Higgs doublet model contributions and 
minimal supersymmetric extension of the SM (MSSM) to the CP asymmetries were discussed in 
refs. [23 1 and [24|, respectively. Recently, CP violation in the polarized b — » d£ + £~ decay has 
been also investigated in the SM [25 1 and also in a general model independent way [26|. 

The aim of this work is to investigate B — > X^ £ + ^~decay with emphasis on CP violation and 
NHB effects in a CP softly broken 2HDM, which is called model IV in the literature E71E51 . 
In model IV, up-type quarks get masses from Yukawa couplings to the one Higgs doublet, and 
down-type quarks and leptons get masses from another Higgs doublet. In such a 2HDM, all the 
parameters in the Higgs potential are real so that it is CP-conserving, but one allows the real and 
imaginary parts of (j)f(j)2 to have different self-couplings so that the phase £, which comes from 
the expectation value of Higgs field, can not be rotated away, which breaks the CP symmetry (for 
details, see ref ||27K In model IV, interaction vertices of the Higgs bosons and the down-type 
quarks and leptons depend on the CP violating phase £ and the ratio tan = V2/V 1 , where v\ and 
V2 are the vacuum expectation values of the first and the second Higgs doublet respectively, and 
they are free parameters in the model. The constraints on tan f3 are usually obtained from B — B, 
K — K mixing, b — > S7 decay width, semileptonic decay b — > ctv and is given by [29| 

0.7<tan/?<0.52(-^-), (1) 
1 GeV 

and the lower bound m H ± > 200 GeV has also been given in [29 1. As for the constraints on £, 
it is given in ref. [ 27 ] that y/\ sin 2£| tan f3 < 50, which can be obtained from the electric dipole 
moments of the neutron and electron. 

For inclusive B-decays into lepton pairs, in addition to the CP asymmetry and the forward- 
backward asymmetry, there is another parameter, namely polarization asymmetry of the final 
lepton, which is likely to play an important role for comparison of theory with experimental 
data. It has been already pointed out 1 30 1 that together with the longitudinal polarization, Pl, 
the other two orthogonal components of polarization, transverse, Pj>, and normal polarizations, 
P/v, are crucial for the t + t~ mode since these three components contain the independent, but 
complementary information because they involve different combinations of Wilson coefficients in 
addition to the fact that they are proportional to me/mb. 
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The paper is organized as follows: Following this brief introduction, in sectional we first 
present the effective Hamiltonian. Then, we introduce the basic formulas of the double and dif- 
ferential decay rates, CP violation asymmetry, Aqp, forward-backward asymmetry, Afb, and 
CP violating asymmetry in forward-backward asymmetry Acp(Afb) for B — ► £ + £~ decay. 
Section|3]is devoted to the numerical analysis and discussion. 



2 The Effective hamiltonian for B -> X d £ + £- 

It is well known that inclusive decay rates of the heavy hadrons can be calculated in the heavy 
quark effective theory (HQET) [31 1 and the important result from this procedure is that the lead- 
ing terms in 1/ m q expansion turn out to be the decay of a free quark, which can be calculated in 
the perturbative QCD. On the other hand, the effective Hamiltonian method provide a powerful 
framework for both the inclusive and the exclusive modes into which the perturbative QCD cor- 
rections to the physical decay amplitude are incorporated in a systematic way. In this approach, 
heavy degrees of freedom, namely t quark and W^,H^, h°, H° bosons in the present case, are 
integrated out. The procedure is to take into account the QCD corrections through matching the 
full theory with the effective low energy one at the high scale n = m\y and evaluating the Wilson 
coefficients from m\v down to the lower scale \i ~ 0(m&). The effective Hamiltonian obtained 
in this way for the process b — > d£ + £~ , is given by fT9ll20l : 

AC ( 10 10 

H eff = -l-WECiMOiM + J^QiM 

V 2 I i=i i=l 

- \ u {C x {n)p^) - OiOz)] + C 2 Qj)[0$(jjl) - 2 ( M )]}| (2) 

where 

y tbV td 

and we have used the unitarity of the CKM matrix i.e., V t bV t * d + V u t,V* d = —V c bV* d . The explicit 
forms of the operators Oj can be found in 0. 0\ and are the new operators for b — > d 
transitions which are absent in the b — ► s decays and given by 

O u x = (d almu P L up){u pl mu P L d a ) 
O2 = {d almu PLU a ){u pl mu PLdp). 

The additional operators Qi(l = 1, .., 10) come from the NHB exchange diagrams and are de- 
fined in ref. lTl9l . 

In Eq.©, Ci(fi) are the Wilson coefficients calculated at a renormalization point fi and their 
evolution from the higher scale \i = mw down to the low-energy scale \x = is described by 
the renormalization group equation. Although this calculation is performed for operators 0{ in the 
next-to-leading order (NLO) the mixing of Oi and Qi in NLO has not been given yet. Therefore 
we use only the LO results. The form of the Wilson coefficients d(mb) and Cg^mi,) in the LO 
are given in refs. HJ and |[T9ll27ll . respectively. 

We here present the expression for Cg(/i) which contains, as well as a perturbative part, a part 
coming from long distance (LD) effects due to conversion of the real cc into lepton pair £ + £~: 

C 9 e// (/u) = C$ ert (v) + Y resm (s) , (4) 
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where 



Cf rt (//) = C 9 + h(u,s)[3C 1 (ri + C 2 (fi) + 3C 3 (ri + C^)+3C 5 ( f i) + C 6 (ri 

+ A„(3d + C 2 )\ - \h{l, s) {ACM + 4C 4 (/u) + 3C 5 (m) + C 6 0)) 

- ^(M [C 3 (» + 3C 4 (aO + A u (6Ci(/i) + 2C 2 (m))] 

+ ^(3C3(/i) + C4(/i)+3C 5 (//) + C 6 (/i)) , 



(5) 



and 



.to 



^ 7rr(Vj -^+r>^ 



a 



x [(3Ci(/x) + C 2 (m) + 3C 3 (^) + C 4 (/z) + 3C 5 (m) + CM) 

+ X U (3CM + C 2 (fx))} . (6) 



In Eq.©, s = q 2 /m 2 B where q is the momentum transfer, 



nib 



and the functions h(u, s) arise 



from one loop contributions of the four-quark operators 0\ — 0% and are given by 



h(u, s) 



8 m b 
■- In — 

9 n 



— In it H 1 — y 

9 27 9 s 



(7) 



-n( 2 + 2/)l 1 -y 



1/2 



In 



n=y±i 



itt), fory = ^<l 



2 arctan • 



for y = ^- > 1 



WO.s = in msH — ?7r 

V ' ; 27 9 n 9 9 



(8) 



The phenomenological parameter k in Eq. © is taken as 2.3 (see e.g. [32 1). 

Next we proceed to calculate the differential branching ratio dBR/ds, forward-backward 
asymmetry Afb, CP violating asymmetry Aqp, CP asymmetry in the forward-backward asym- 
metry Acp{Afb) and finally the lepton polarization asymmetries of the B — » £ + ^~decays. 
In order to find these physically measurable quantities we first need to calculate the matrix ele- 
ment of the B — > X^l + i~ decay. Neglecting the mass of the d quark, the effective short distance 
Hamiltonian in Eq.© leads to the following QCD corrected matrix element: 



M 



GpOi 



VtbV t U C e 9 ff d 7M (l - 75)&#y^ + Cio d 7M (l " K)bljV 



- 2C e / f ^dia^q'il + lB )bhH + C Ql d(l + l5 )bU + Cq 2 d(l + 75 )&«| -(9) 

When the initial and final state polarizations are not measured, we must average over the initial 
spins and sum over the final ones, that leads to the following double differential decay rate 



3a 2 



\V tb V, 



* |2 
td\ 



d 2 T 

Ud~z = r(jB ?_> ' Xce ^'47r 2 f(u)k(u)^ " J \V cb \ 2 

+ 2 (l + ^) Re(C e 7 ff C e 9 ff *) + v S zRe(C 10 C e 9 ff *) 

+ ^Re((2C 7 e// + C e 9 ff )C* Ql ) + VtRc(C w C* Q2 ) 



v[2vzRt{C e 7 fS C* lQ ) 



+ 



+ 



(l + s)-{l-s)v 2 z 2 + 4t 



\C 



e//|2 



At 

s 



\C 



e//|2 



(1 + a) - (1 - s) v 2 z 2 - it] |do| 2 + -s\C Q2 \ 2 + -{s- At)\C Ql \ 2 } (10) 



where v = y/1 — 4i/s, t = mf/mf and z = cos 6, where 9 is the angle between the momentum 
of the B-meson and that of t~ in the center of mass frame of the dileptons £~£ + . In Eq. dTOl . 

T(B -+ X c lv) = ( §^t\V cb \ 2 f(u)k(u) , (11) 

where 

f(u) = 1 - 8u + 8u 4 - u 8 - 24u 4 ln(u) (12) 

(13) 



k(u) = 1 - 2a ° {mb) 



3tt 



2 31\ , n a2 , 3 
vr 2 - T J(l-m c 2 ) + - 



are the phase space factor and the QCD corrections to the semi-leptonic decay rate, respectively, 
which is used to normalize the decay rate of B — > X d £ + £~ to remove the uncertainties in the 
value of nib- 

Having established the double differential decay rates, let us now consider the forward-backward 
asymmetry Ap b of the lepton pair, which is defined as 



a i \ JO dsdz J — l dsdz ci /i\ 

Afb ^ = ri, d 2 T , r0 T- d*r • ( 14 ) 



JO uz dsdz J -I az dsdz 

[ l dz^- + f° dz-^- 

JO " Z (isd2 ^ J-l " Z dsd2 

The Af\b's for the B — > X d £ + ^~decays are calculated to be 



A FB (s) = J^Re[C w (2C e 7 ff + sC e 9 ff *)] + VtRe[C Ql (2C 7 e// * + C e 9 ff *)], (15) 

where 

A( S ) = (S + + 2t ~ 8st) |C 10 | 2 + - 2 {2 + s)(s + 2t)|C 7 ^| 2 + (1 + 2s)(l + -)|C 9 e// | 2 
s s^ s 

+ H( s + 2t)Re(C 7 e// C 9 e// *) + 6^Re(C 9 e// ^ 2 ) + |(a - 4i)|C Ql | 2 + | S |Cq 2 | 2 , 

(16) 

which agrees with the result given by ref. |U, in case of switching off the NHB contributions and 
setting rri£ = 0, but differs slightly from the results of l24l . 

We next consider the CP asymmetry Acp between the B — > X<i& + £~ and the conjugated one 
73 — ► , which is defined as 

dT _ d£ 

= frf < 17 > 

ds ds 

where 

dr _ dr{B -» dr _ dr(ff^x d i+r) 

' ds ds 

After integrating the double differential decay rate in Eq.dTDTl over the angle variable, we find 
for the B X d £ + £~ decays 

f-^-^'j^ky^-^w^ 1 ?^'- (19) 

For the antiparticle channel, we have 



We have also a CP violating asymmetry in A F b, Acp(Apb), in B — > £ + £~ decay. Since 
in the limit of CP conservation, one expects Aps = —Afb EEJ, where App and App are the 
forward-backward asymmetries in the particle and antiparticle channels, respectively, Acp(App) 
is defined as 



A c p (A FB ) 



Apn + A 



FB 



with 



Afb = A FB {\ U — > A*;£ — * — £) 



(21) 



(22) 



Finally, we would like to discuss the lepton polarization effects for the B — > X& £ + £ decays. 
The polarization asymmetries of the final lepton is defined as 



Pjs) 



(dT(S n )/ds) - (dT(-S n )/ds) 
(dT(S n )/ds) + (dT(-S n )/ds) ' 



(23) 



for n = L, N, T. Here, Pl, Pt and P/v are the longitudinal, transversal and normal polarizations, 
respectively. The unit vectors S n are defined as follows: 



Sl 
Sn 
St 



(o,&) : 

(0,ejv) 
(0,e T ) 



0, 



P- 







IP- 1 / 
p x p_ 
\p x P- 1 

0, 6jv x e L 



(24) 



where p and p_ are the three-momenta of d quark and ^ lepton, respectively. The longitudinal 
unit vector Sl is boosted to the CM frame of £ + £~ by Lorentz transformation: 



Sl,cm 



\P- 1 E t p- 



m e ' mi\p-\ 



(25) 



It follows from the definition of unit vectors S n that Pt lies in the decay plane while Pn is 
perpendicular to it, and they are not changed by the boost. 

After some algebra, we obtain the following expressions for the polarization components of 
the £~ lepton in B —> X^ £ + £~ decays: 



Pl 
Pt 



A 



2Ci (6C 7 e// '* + (l + 2s)C 9 
3ttV7 / 4 



3C Ql (2 Vt C w + sC* Q2 



C e 7 ff \ 2 -s\C e 9 ff \ 2 +Re 



2C 7 e// *(Ci -2(7 9 e// * 



+ 



2-Jl 



C5J 



C '9 // (ClO + -^^ 2 ) + 



■4t 



2\/i 



2y/i 



(26) 



3irv 



Im 



C 10 ( S C^ 2 + 2V~t(C e / f * + sC e 9 ff *)) + sC Ql (2C e / f * + CI"*] 
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3 Numerical results and discussion 



In this section we present the numerical analysis of the inclusive decays B — > X^i + l'm model 
IV. We will give the results for only I = r channel, which demonstrates the NHB effects more 
manifestly. The input parameters we used in this analysis are as follows: 

m b = 4.8 GeV , m c = 1.4 GeV , m t = 175 GeV , m T = 1.78 GeV , 

BR{B -> X c ev e ) = 10.4% , m H ± = 200 GeV , m H o = 160 GeV , m h o = 115 GeV 

a' 1 = 1/129 , G F = 1.17 x 10~ 5 GeV~ 2 . (27) 

The Wolfenstein parametrization [ 34 1 of the CKM factor in Eq. © is given by 

A„= "VV! 2 V" +0(A 2 ), (28) 

(i - pY + if 



and also 



\v tb v t 



* |2 



A 2 [(l-p) 2 + 7 ? 2 ] + 0(A 4 ). (29) 



The updated fitted values for the parameters p and i] are given as l35l 

p = 0.22 ±0.07 (0.25 ±0.07), 

fj = 0.34 ± 0.04 (0.34 ± 0.04) , (30) 

with (without) including the chiral logarithms uncertainties. In our numerical analysis, we have 
used (p, rj) = (0.25; 0.34) . 

The masses of the chained and neutral Higgs bosons, m H ±, m H o, and m h o, and the ratio 
of the vacuum expectation values of the two Higgs doublets, tan 0, remain as free parameters 
of the model. The restrictions on m#±, and tan/? have been already discussed in section [2 
For the masses of the neutral Higgs bosons, the lower limits are given asm^o > 115 GeV and 
m h o > 89.9 GeV in ESI. 

In the following, we give results of our calculations about the dependencies of the differ- 
ential branching ratio dBR/ds, forward-backward asymmetry Afb(s), CP violating asymme- 
try Acp(s), CP asymmetry in the forward-backward asymmetry Acp(Afb){s) and finally the 
components of the lepton polarization asymmetries, Pl(s), Pt(s) and P^(s), of the B — > 

r + r~decays on the invariant dilepton mass s. In order to investigate the dependencies of 
the above physical quantities on the model parameters, namely CP violating phase £ and tan (5, 
we eliminate the other parameter s by performing the s integrations over the allowed kinematical 
region so as to obtain their averaged values, < Afb >. < Acp >, < Acp{Afb) >, < Pl >, 
< P T > and < P N >. 

Numerical results are shown in Figs. (TTb- dT3T > and we have the following line conventions: 
dashed lines, dot lines and dashed-dot lines represent the model IV contributions with tan = 
10, 40, 50, respectively and the solid lines are for the SM predictions. The cases of switching off 
NHB contributions i.e., setting Cq z = 0, almost coincide with the cases of 2HDM contributions 
with tan /? = 10, therefore we did not plot them separately. 

In Fig.Q, we give the dependence of the dBR/ds on s. From this figure NHB effects are 
very obviously seen, especially in the moderate-s region. 

In Fig. (0 and Fig. (J3J, Afb(s) and < Afb > as a function of s and CP violating phase £ 
are presented, respectively. We see that Afb is more sensitive to tan (3 than the dBR/ds and it 
changes sign with the different choices of this parameter. It is seen from Fig.© that < Afb > 
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is quite sensitive to £ and between (0.15,0.28) x 10 _1 . We also observe that < App- > differs 
essentially from the one predicted by the CP-conservative 2HDM (model II, for examples, see 
[37 1), which is 0.028 and 0.023 for tan/? = 40, 50, respectively. In region 1 < £ < 2 change in 

< Afb > with respect to model II reaches 25%. 

Figs. © and ^ show the dependence of Acp(s) on s and < Acp > on £, respectively. We 
see that Acp{s) is also sensitive to tan f3 and its sign does not change in the allowed values of s 
except in the resonance mass region. It follows from Fig. (|5} that < Acp > is not as sensitive as 

< App > to £, and it varies in the range (0.15, 0.33) x 10 _1 . 

Acp(Afb)(s) and < Acp{App) > of B — > XdT + r~as a function of s and CP violating 
phase £ are presented in Fig. © and Fig. (Q, respectively. We see that Acp(App)(s) changes 
sign with the different choices of tan/3. < Acp{App) > is between (0.010,0.040) and differs 
essentially from the one predicted by model II, which is 0.038 and 0.027 for tan/? = 40,50, 
respectively. In region 1.5 < £ < 2.5 change in < App > with respect to model II reaches 35%. 

In Figs. (l8l- (fT0t . we present the s dependence of the longitudinal Pp, transverse Pp and 
normal P/v polarizations of the final lepton for B — » Xd t + t~ decay. It is seen that NHB con- 
tributions changes the polarization significantly, especially when tan (3 is large. We also observe 
that except the resonance region, Pp is negative for all values of s, but Pp and P/v change sign 
with the different choices of the values of tan (3. In Figs. (fm> - (fT3l . dependence of the averaged 
values of the longitudinal < Pp >, transverse < Pt > and normal < P/v > polarizations of 
the final lepton for B — ► X^r + r^ decay on £ are shown. It is obvious from these figures that 

< P/v > and < Pp > are more sensitive to £ than < Pp >. In region 1.5 < £ < 2.0 change in 

< P/v > with respect to model II reaches 25%. Thus, measurement of this component in future 
experiments may provide information about the model IV parameters. 

Therefore, the experimental investigation of App, Acp, Acp{App) and the polarization 
components in B — > Xd decays may be quite suitable for testing the new physics effects 
beyond the SM. 
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Figure 3: < A FB > as a function of £. 
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Figure 4: The CP asymmetry as a function of s, where £ = tt/4. 
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Figure 6: The CP asymmetry in the forward-backward asymmetry as a function of s, where £ = 7r/4. 
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Figure 7: The CP asymmetry in the forward-backward asymmetry as a function of £. 
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Figure 9: Pt(s) as a function of s , where £ = tt/4. 
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Figure 10: Pn(s) as a function of s , where £ = 7r/4. 
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Figure 1 1: < P L > as a function of £. 
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Figure 12: < P T > as a function of £. 
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Figure 13: < P N > as a function of £. . 
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